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We study the three dimensional SU(2)-symmetric noncompact CP^ model, with two charged 
matter fields coupled minimally to a noncompact Abelian gauge-field. The phase diagram and the 
nature of the phase transitions in this model have attracted much interest after it was proposed to 
describe an unusual continuous transition associated with deconfinement of spinous. Previously, it 
has been demonstrated for various two-component gauge theories that weakly first-order transitions 
may appear as continuous ones of a new universality class in simulations of relatively large, but 
finite systems. We have performed Monte-Carlo calculations on substantially larger systems sizes 
than those in previous works. We find that in some area of the phase diagram where at finite 
sizes one gets signatures consistent with a single first-order transition, in fact there is a sequence 
of two phase transitions with an 0(3) paired phase sandwiched in between. We report (i) a new 
estimate for the location of a bicritical point and (ii) the first resolution of bimodal distributions 
in energy histograms at relatively low coupling strengths. We perform a flowgram analysis of the 
direct transition line with rescaling of the linear system size in order to obtain a data collapse. The 
data collapses up to coupling constants where we find bimodal distributions in energy histograms. 
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I. INTRODUCTION 

Recently, the CP^ model consisting of two matter 
fields coupled to an Abelian gauge field has been of 
great interest in condensed matter physics. One of the 
sources of interest is the proposed concept of deconfined 
quantum criticality (DQC). It has been intensively de- 
bated as a possible novel paradigm for quantum phase 
transitions .'^"i^ Such quantum criticality has been sug- 
gested to describe phase transitions that would not fit 
into the Landau-Ginzburg- Wilson (LGW) paradigm of a 
continuous (second-order) phase transition.""'^ In partic- 
ular, the continuous quantum phase transition from an 
antiferromagnetic Neel st ate in to a paramagnetic valence- 
bond solid (VBS) state|--'-°' does not agree with the 
LGW-description, according to which two phases with 
different broken symmetries generically are separated by 
a first-order phase transition. Recently, evidence for the 
DQC scenario has been claimed in studies of the so-called 
J-Q model, ^ which is a Heisenberg model with additional 
higher-order spin interaction terms. Namely, it was sug- 
gested that high-precision Quantum Monte Carlo sim- 
ulations of this model support a continuous Neel - VBS 
phase transition in accordance with the DQC scenario.'^HIl 

It has been proposed that the critical field theory of a 
continuous Neel - VBS phase transition is the so-called 
noncompact CP^ model (NCCP^), with a SU(2) sym- 
metric field coupled to a no ncompact U(l) gauge field in 
three dimensions (3D) Initial efforts on studying this 
effective model were focused on the special case where 
the SU(2) symmetry was broken down to a U(l)xU(l) 
symmetry, i.e., the easy-plane limit. For this case, a 
continuous phase transition was claimed.™ However, in 
Ref.lTOl the existence of a paired phase in the U(l)xU(l) 
easy-plane action was pointed out. (For earlier discus- 



sions of paired phases in various U(l)xU(l) systems, see 
Refs. [TTJ [HIH2S1) Furthermore, resorting to mean-field 
theory arguments, it has been pointed out that at least 
in the vicinity of a paired state (in the parameter space of 
the model), the direct phase transition from a symmet- 
ric state to a state with broken U(l)xU(l) symmetry, 
should be first-order.-" Subsequent Monte-Carlo calcu- 
lations have reported a weak firs t-ord er phase transition 
for the easy-plane NCCP^ vaodelW^The so-called fiow- 
gram method has also been introduced in Ref. [TOl specif- 
ically to characterize weak first-order phase transitions. 
Using this method, the direct phase transition from a 
symmetric state to a state with broken U(l)xU(l) sym- 
metry, has been claimed to be first-order for any non-zero 
value of the coupling constant. The phase transitions in 
the easy-plane limit of the NCCP^ model were also exten- 
sively studied in variety of other regimes in the context 
of two-component su percon ductors with independently 
conserved condensates ! ^ ^ 1 ^ ^I MI 

For the SU(2)-symmetric case, Monte Carlo compu- 
tations have been performed in Ref. 9. Here, a direct 
second-order phase transition was suggested, but the sys- 
tem sizes that were considered were quite small. In a 
subsequent paper,'^^! an extensive study of the model was 
performed. In particular, for the direct transition line, 
a second-order phase transition was claimed. At higher 
couplings to the gauge field, it was suggested to turn 
into a first-order transition via a tricritical point. On the 
other hand, in Ref. (see also Ref. "Tf), it was argued 
that the direct transition line is first-order. The flowgram 
method employed in Ref. 1 141 showed no evidence for a tri- 
critical point along the direct transition line. Rather, in 
this work the large-scale behavior at small couplings to 
the gauge field was found to be the same as for higher cou- 
plings, where indications of a first-order transition were 
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seen by resolving a bimodal distribution in the energy his- 
tograms. Note that at the system sizes studied in Ref.fT4l 
no bimodal distributions were resolved at small coupling 
constants. Nonetheless, in Rcf. TT, it was concluded that 
even for weaker couplings, bimodal distributions indica- 
tive of a first-order phase transition would emerge for 
large enough system sizes. This conclusion was based on 
the similarity of scaling of various quantities for large and 
small couplings, as evidenced by the flowgrams. 

The weakness of the observed first-order phase transi- 
tion, combined with the necessity of assessing the order 
of the phase transitions also in the limit of vanishingly 
small coupling strength, renders this problem computa- 
tionally extremely demanding. In this work, we therefore 
examine the phase diagram of the NCCP^ model at sub- 
stantially larger systems sizes than what has been done 
in previous worksI^^Eil Performing the computations on 
larger systems allows us to perform a very detailed inves- 
tigation of the range of parameters where a paired phase 
is sandwiched in between the fully disordered and fully 
ordered state. This means that these two phases are sep- 
arated, not by a direct transition, but by two separate 
transitions. At small system sizes, these two separate 
phase transitions in fact give signatures which would lead 
one to conclude that the system features one single first- 
order phase transition. The existence of a paired phase 
sandwiched in between the fully ordered and disordered 
states emerges only when one considers large enough sys- 
tems. Our study thus allows us to provide improved es- 
timates for the location of a bicritical point where the 
direct transition line splits into two. 



II. MODEL 

The continuum NCCP"'^ model is written as 



H 



(1) 

y - ieA(x)] *(x)|' -t- [V X A(x)]'} , 

(2) 



where /3 is the inverse temperature and 'J'^(x) = 
(V'i(x), V'2(x)) are two complex fields that are coupled to 
a noncompact gauge field A(x) with charge e. The fields 
■0c(x), c e {1,2}, obey the CP^ constraint, |*(x)| = 1. 

The model can be mapped onto a nonlinear 0(3) a 
model coupled to massive vector fields.l^Sl By introducing 
the fields. 



C(x) ^'-Y. [V'c(x)VV':(x) - ^:(x)VVc(x)] - eA(x) 



n(x) = \I'"^(x)cr*(x), 



(3) 
(4) 



where the components of cr are the Pauli matrices, the 
NCCP^ model (fll) can be rewritten aJ^l 



i/=i[a,n(xr + i[C(xf 




9.Ca(x) 



n(x) • 5^n(x) x d\n{y.) 



(5) 



where sum over repeated indices is assumed. The model 
represents an 0(3) nonlinear a model coupled to a mas- 
sive vector field C(x). The latter represents a charged 
mode, and its mass is the inverse magnetic field pene- 
tration length. At least for sufficiently large values of 
electric charge coupling, the model can undergo a Higgs 
transition (where gauge field becomes massless) without 
restoring simultaneously any broken global symmetries. 
In that case, the remaining broken global symmetry is 
0(3) which is described by the order parameter n(x). 

If one introduces an easy-plane anisotropy for the vec- 
tor field n(x), this would break the symmetry of the 
model to U(l)xU(l), and the separation of variables 
yields a neutral and a charge d rnode, the physics of which 
has been extensively studied.fiSinillHll! However, there is 
one substantial difference in the case of SU(2) symmetry. 
The charged and neutral sectors are coupled through the 
last term in Eq. ([5|. Another difference compared to 
the U(l)xU(l) case is that in two dimensions, stable 
singly quantized vortex lines do not exist in a type-II 
SU(2) model (the same applies to vortex lines in three 
dimensions). l^S! On the other hand, a type-I SU(2) model 
has energetically stable counterparts of ordinary singly 
quantized type-I vortices. Since composite vortices are 
topological excitations which lead to the occurrence of 
paired states in U(l)xU(l) systems, this aspect makes 
the phase diagram of SU(2) theory an especially inter- 
esting problem to study. 

In the Monte Carlo simulations, we employ a lattice 
realization of this model on a cubic lattice with size 
and with lattice constant a — 1. The fields ■(/'c(x) are 
then defined on the vertices r e {ix -I- jy -I- kz\i,j, k g 
{1, . . . , L}} of the lattice, '0c(x) — ipc.r- For the first 
term in (pi), we rescale the gauge field by and invoke 
the gauge invariant lattice difference. 



d 
dx,, 



ie^p(x) 



V'c(x) -J> ip. 



(6) 



where fi G {x, ?/, z} and r + p, denotes the nearest- 
neighbor lattice point to vertex r in the /i-direction. The 
gauge field A^^ ,. lives on the {r,r + p,) links of the lattice. 
For the Maxwell term we get 

[V X A(x)]^ ^ J2 e^-AA.AA.r, (7) 



where is the forward finite difference operator, 
AjyAxr = A\^r+o — A\^r, and e^i^A is the Levi-Civita 
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symbol. In addition, by invoking the CP^ constraint and 
discarding constant factors in the partition function Z, 
we obtain the following lattice realization of the NCCP^ 
model: 

//>1 /•27T />27r 

VA Vu VOi / Ve-2 e-^", (8) 
Jo Jo Jo 



- VUr^/Ur+p. cos (A^^i^r - A^^r) 

- Vl - UrV^l - "r+A COS {A^02,r " ^p,r 



2e2 



where Ur = IV'i.rl^ = 1 ~ |V'2.r|^ and where |V'c,r| is the 
amplitude and 0c,r is the phase of the complex fields ipc.r- 



III. DETAILS OF THE MONTE CARLO 
SIMULATIONS 

The Monte Carlo simulations are performed on a cubic 
lattice with periodic boundary conditions in all directions 
and with size where L G {8, . . . , 96}. Up to 4.0 • 10^ 
sweeps over the lattice were performed for the largest 
systems, while up to 1.0 • 10^ sweeps were used for initial 
equilibration and initialization of the coupling distribu- 
tion (see below). Monte-Carlo time-series were routinely 
inspected for equilibration. To test for ergodicity, typi- 
cally 4 independent large simulations were performed for 
the largest system sizes. Histograms based on raw data 
and reweighted data were also compared for consistency. 
For most of the simulations, the parallel tempering (PT) 
algorithm was employed.l^HHSi] To be specific, we fix the 
coupling e and perform the computations on a number of 
replicas (typically from 8 to 32 depending on the system 
size L and the range of /3 values) in parallel at differ- 
ent values of /3. A Monte Carlo sweep consists of sys- 
tematically traversing all lattice points with local trial 
moves of all six field variables by the Metropolis-Hastings 
algorithm.i^^EII For u^, the proposed new values are cho- 
sen with uniform probability within the interval [0,1], 
and for 0c,r, the proposed new values are chosen with 
uniform probability within the interval [0,27t). For the 
noncompact gauge field, the proposed new values are cho- 
sen within some limited increment (typically [— 7t/4, 7t/4]) 
from the old values.'^ There is no gauge fixing involved 
in the simulations. In addition to these local trial moves, 
the Monte Carlo sweep also includes a PT trial move of 
swapping replicas at neighboring /3 values. 

All replicas were initially thermalized from an ordered 
or disordered start configuration. Then, initial runs were 
performed in order to produce an optimal distribution of 
couplings for the simulation. In some cases, the set of 
couplings was found by measuring first-passage-times.'^ 
In this approach, the optimal set of couplings maximizes 



the flow of replicas in parameter space, essentially by 
shifting coupling values towards the bottlenecks.'^ How- 
ever, in cases with no severe bottleneck, the optimal set 
of couplings was found by demanding that the acceptance 
rates for swapping neighboring replicas were equal for all 
couplings.'^ Irrespective of how the set of couplings was 
found, it was always ascertained that replicas were able to 
traverse parameter space sufficiently many times during 
production runs. The measurements were postprocessed 
by multiple histogram reweighting.^Ss! Random numbers 
were generated by the Mcrscnnc-Twister algorithm.!^ Er- 
rors were determined by the jackknife method.'^ 

As mentioned in the introduction, the NCCP^ model is 
a difficult model on which to perform Monte Carlo com- 
putations. In Ref. T?', the NCCP^ model was mapped 
to a so-called J-current mode l, wh ich allows simulations 
based on the worm algorithm.SlEIl (por the sake of com- 
pleteness, and since to our knowledge the details of the 
mapping have not been published, we present the deriva- 
tion of this mapping in Appendix|A]) An approach based 
on the J-current model was attempted as well. However, 
due to the presence of long-range interactions in this for- 
mulation, it was difficult to work with the lattice sizes 
above L ^ 40. Hence, the computations were performed 
on the model in the original NCCP^ formulation, using 
the PT algorithm with which it is easy to grid-parallelize 
the lattice. 



IV. OBSERVABLES AND FINITE-SIZE 
SCALING 

Perhaps the most familiar quantity that is used to ex- 
plore phase transitions, is the specific heat Cy. The spe- 
cific heat is given by the second moment of the action. 



(9) 



where brackets (...) denote statistical averages. In most 
cases, Cy exhibits a well-defined peak at the phase tran- 
sition. For a continuous phase transition the correlation 
length diverges with critical exponent ly as ^ ^ \t\ , 
with t — {j3 — /3c) //3 being the deviation from the critical 
coupling /3c. The critical exponent a is defined by the 
singular part of Cy, given by Ct, ^ |i| . Then, in a 
limited system of size , the finite-size scaling (FSS) of 
the specific heat is given by 



Cy^C0 + Cli"/^ 



(10) 



where Co and Ci are non-universal coefficients. For 
a first-order transition, with two coexisting phases and 
no diverging correlation length, there is indeed no crit- 
ical behavior. Still, first-order transitions exhibit well- 
behaved FSS with "effective" exponents, a = 1 and 
1/ = 1/3.1321331 Hence, the peak of the specific heat scales 
as 



Cy 



(11) 
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for a first-order transition. Distinguishing between con- 
tinuous and first-order transitions is an important issue 
in tlie present work. For that purpose, FSS of the specific 
heat peak will play an important role. 

We a lso in vestigate the third moment of the action 

given bj^mni 



M,^^{{H-{H)y 



(12) 



In the vicinity of the critical point, this quantity typically 
features a minimum point and a maximum point [see for 
instance the inset in panel (b) of Fig. [2] . The difference 
in the M3 value of these two extrema scales as 



(AM; 



height 



(13) 



and the difference in the coupling values scales as 

(AAf3)width - i-'/", (14) 

for a continuous phase transition. For a first-order tran- 
sition, the FSS is 



(AM3)height 



and 



(AM3), 



idth 



(15) 



(16) 



As was mentioned above, one may construct a three 
component gauge neutral field rir ,[231 given by 



rir = *;0-*r, 



(17) 



where the components of cr are the Pauli matrices. Since 
it is a unit 0(3) vector, we can introduce a "magnetiza- 
tion" , 



M = ^rir 



(18) 



The order parameter (m), where m = M/L^, signals 
the onset of order in the 0(3) gauge neutral vector field 
rir, and the critical point of this transition can be accu- 
rately determined by a proper analysis of the finite-size 
crossings of the associated Binder cumulantjI^^SH 



3(M4) 
2(M2)2- 



(19) 



The finite-size crossings of the Binder cumulant are 
known to converge rapidly towards the critical coupling 
/3c. Hence, /3c can be accurately determined by a simple 
extrapolation of the finite-size crossings to the thermo- 
dynamic limit or by invo king scaling forms that account 
for finite-size corrections .H2Eni 

A number of quantities related to magnetization may 
be used to extract critical exponents from the Monte 
Carlo simulations. The magnetic susceptibility, given by 



X = L3/?(m2), 



(20) 



when /3 < /?c, scales as x '-^ L^-i) p — Hence, 
we may determine the anomalous scaling dimension 77 by 
FSS of X measurements obtained at /3c- 

The exponent v can, alternatively, be determined by 
calculating the logarithmic derivative of the second power 
of the magnetization,'^ 



d_ 
d(3 



In (to^) 



(H). 



(21) 



The FSS of this quantity is ^ln{m'^) - L^/'^. Since 
the logarithmic derivative exhibits a peak that is associ- 
ated with the critical point, it is possible to extract v by 
measuring the logarithmic derivative at the pseudocriti- 
cal point, without an accurate determination of /3c. 

Similar to Ref. [T31 we search for the critical point of 
the Higgs transition by measuring the dual stiffness 



PSuai(q) 



Er,.,AeA.^AA.A;,,re"l"- 



(271)2^3 



(22) 



which is the Fourier space correlator of the magnetic 
field. This order parameter for the Higgs transition is 
dual in the sense that it is finite in the high-temperature 
phase and zero in the low-temperature phase. Like in 
Ref. 1131 this quantity is measured at the smallest avail- 
able wavevector q ^ 0. We chose to measure p^uai 
Qmin = (27t/L,0,0). At the critical point, the quantity 
^PduaiC^min) Universal, such that the finite-size cross- 
ings of ip[Juai('^min) '^^^ tie uscd to estimate the critical 
point of the Higgs transition. In addition, measuring the 
coupling derivative of Lp'^'^^^{qj-^^jJ can be used to esti- 
mate the correlation length exponent v as 



_d_ 

a/3 

at the critical point. 



(23) 



V. NUMERICAL RESULTS 

A. Outline of the phase diagram 

The phase diagram of the NCCP^ model is presented 
in Fig. [1] For small values of /3, there is a normal phase 
that can be recognized by a disordered gauge neutral vec- 
tor field rir and a massless gauge field. Hence, (to) — 
and Pduai(l) 7^ in this phase. For large values of e 
and higher values of /3, there is a transition into a phase 
that we label the 0(3) phase. Here, the vector field 
rir is ordered (the 0(3) symmetry is spontaneously bro- 
ken), (to) ^ 0, whereas the gauge field remains massless, 
/^duai(l) 7^ 0. In the case of U(l) x U(l) symmetric su- 
perconductors, this phase is sometimes denoted a metal- 
lic superfluid or a paired phase, with long-range order 
in the gauge neutral linear combination of the phases 
(in the U(l) x U(l) case), but not in the individual 
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Figure 1. (Color online) Phase diagram of the NCCP^ model. 
SU(2) phase: Fully ordered phase where the 0(3) symme- 
try is spontaneously broken, (m) 7^ 0, and the gauge field is 
massive, Pduai(l) = 0. 0(3) phase: 0(3) symmetry is spon- 
taneously broken, (m) ^ 0, but the gauge field is massless, 
Pduai(l) 7^ 0. Normal phase: 0(3) symmetry is restored, 
(m) = 0, and the gauge field is massless, p^uai(l) 7^ 0- '^^^ 
rect transition line from the SU(2) phase to the normal phase 
is denoted by -f- markers and a solid red line. The Higgs 
transition line between the SU(2) phase and the 0(3) phase 
is denoted by *-markers and a dotted blue line. The transition 
line between the 0(3) phase and the normal phase is denoted 
by X -markers and a dashed green line. Lines are guide to the 
eyes. 



o^eg| io | ii | 2iH23 | 26 | 52 | p^.^^^ ^j^g o(3) phase, by reducing 
the value of e, one enters an ordered phase that we label 
the SU(2) phase. Going into this phase, the gauge field 
dynamically acquires a Higgs mass and the system be- 
comes a two-component NCCP^ superconductor. Note 
that the Higgs transition is related to a local symmetry, 
and indeed, is not associated with spontaneous symme- 
try breaking.'^ This aspect should be kept in mind where 
we for brevity refer to the fully ordered state as "broken 
SU(2)" or "fully broken state" to distinguish it from a 
paired state. The SU(2) phase is recognized by measur- 
ing (m) ^ and p^ZJ^^) = 0. 



It is generally expected that at small values of e, the 
SU(2) phase may also be entered directly from the nor- 
mal phase, i.e., without going through the intermediate 
paired phase. The nature of the phase transition along 
this direct transition line in this and related multicom- 
ponent models has been intensively debated due to its 
relevance to deconfined quantum criticality. We will re- 
turn to the direct transition line in Sees. IV Bl and IV CI 
First, we present results for the two separate transition 
lines. 



1. 0(3) line 

In Refs. [T31 EH and [S5J the existence of an inter- 
mediate paired phase, separating a fully ordered state 
from a fully disordered one, was shown in the SU(2)- 
symmetric theory. The nonlinear a model mapping pre- 
sented above, suggests that the transition line between 
the normal phase and the 0(3) phase should be a contin- 
uous transition in the 0(3) universality class, at least in 
the limit far from the bicritical point. We have considered 
this for the case e = 6.0, and the FSS results are given in 
Fig. [2] A log-log plot of the FSS of the peak height in 
In^TO^^ is given in panel (a), and the measured peak 
heights fall on a straight line for L > 20. The best fit to 
the form ^ In (m^) - L^/-" yields 1^ = 0.715 ± 0.004. In 
panel (b), we also measure (AAf3)hoight, and this quan- 
tity exhibits negligible finite-size corrections to scaling 
at least for L > 10. The best fit according to Eq. (13) 
yields a = —0.117 ± 0.011, where the value of 1/ obtained 
above was used. In this case, it was found that v was 
most precisely determined by measuring the peak height 
in J^ln(m^) rather than measuring (AM3)width- The 
maximum peak in M3 is not very sharp [see the inset of 
panel (b)]. Thus, the error bars in (AM3)width are large. 
In order to determine ry, the FSS of the magnetic sus- 
ceptibility X is given in panel (c). Here, x is measured 
at the critical coupling /3c = 2.7894 ± 0.0003, which was 
determined by fitting the Binder crossings of L and L/2 
to a function that accounts for power-law finite-size cor- 
rections. The best fit of x(L) was determined for sizes 
L e {12, . . . , 64} to yield rj = 0.024 ± 0.014. All the ex- 
ponents listed above correspo nd w ell with the exponents 
of the 0(3) universality class .l^^^^ 



2. Superconducting transition. 

Computations have also been performed along the 
transition line between the 0(3) phase and the SU(2) 
phase. In analo gy with the paired phase of the 
U(l)xU(l) modeliojEmHllIIIIll] metallic super- 

fiuid), the transition to the 0(3) sector should be associ- 
ated with the proliferation of single-quanta vortices. In 
the U(l)xU(l) model, such vortices have similar phase 
windings in both complex fields and are topologically 
well-defined objects. In the SU(2) case, such vortices 
can have either similar phase windings in both compo- 
nents, or a phase winding only in one component if the 
other component exists only in the vortex core of the for- 
mer. Such objects are non-topological, and are unstable 
in type-II SU(2) superconductors.^^ This suggests that 
the system should be a type-I SU(2) superconductor in 
order to feature a phase transition into a paired phase. In 
analogy with single-component type-I supercondu ctors, 
one would then expect a first-order phase transition.'^Slini 
A different viewpoint is based on mean-field arguments, 
which suggest that the transition line could be a first- 
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Figure 2. FSS results for the transition between the normal 
phase and the 0(3) phase when e = 6.0. 13 system sizes 
L £ {8, . . . , 64} are used. In all panels, the solid straight 
line is the best fit obtained for a fitting function on the form 
aL*' with two free parameters a and h. Panel (a): Log- 
log plot of the maximum in the logarithmic derivative of 
the second power of the magnetization {d / dl3\n{m?)) p^p.^^ 
[see Eq. (21 1] as a function of L. The best fit is obtained 
for sizes L G {20, ...,64}. The inset shows the measure 
of (9/9/3 ln(m2));3=,3p, in the case when L = 40. Panel 
(b): Log- log plot of the third moment height difference 
(AM3)hcight as a function of L. The best fit is obtained 
for sizes L G {10, ...,64}. The inset shows the measure 
(AM3)hcight in the case when L — 14. Panel (c): Log-log 
plot of the magnetic susceptibility measured at the critical 
coupling X/3=/3c as a function of L. The best fit is obtained 
for sizes L £ {12, . . . , 64}. The inset shows Xi3=)3c f^ir the case 
when L = 40, and the arrowheads indicate that x is measured 
at the same fixed coupling /3c for all sizes. 
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Figure 3. Log-log plot of ^I/p|uai(qmin) nieasured at the 
critical point /3c, as a function of system size L. The charge is 
e = 5.0. Measurements are performed for 15 different system 
sizes L G {8, ...,64}. The derivative was found by calcu- 
lating the differences of Pduai(<lmin)- The solid straight line 
is the best fit obtained with a fit function on the form aL^ 
where a and b are two free parameters. The inset shows the 
^Pduai(lmin) crosslngs for systems L and L/2 as a function of 
L-^. These crossings were used to estimate the critical point, 
/3c = 2.7347 ± 0.0005. Errors in determining /3c are taken into 
account by also considering the sensitivity of v with respect 
to P when estimating the uncertainty in the exponent. 



order transition line in the vicinity of a bicritical point .^^^ 
Other o bjects which can disorder the Higgs sector, are 
Hopfiong^SSni this work we have made no serious at- 
tempts at resolving such topological defects. 

To check the universality class of this line, FSS results 
of 9/9/3 [LpJ^jj[(qj^ijj)], obtained at the critical point with 
e = 5.0, are given in Fig. [3j First, the critical coupling 
was determined to be /3c = 2.7347 ± 0.0005, by consider- 
ing the crossings of L p^^^^y{(i^^^) (see the inset of Fig.|3|. 
Then, the correlation length exponent was estimated to 
he V = 0.664 ± 0.039. This value is consistent with an 
inverted 3Dxy transition line.^^ We have not been able 
to resolve a first-order phase transition at this line. 



B. Estimate for a bicritical point 

In Ref. [131 the flowgram method has been suggested 
as a useful tool to assess whether or not there is a tri- 
critical point at weak couplings to the gauge field. This 
method relies on resolving a first-order phase transition 
at stronger couplings, just below the bicritical point at 
which the paired phase opens up between the normal 
phase and the SU(2) phase. It is thus important to be 
able to determine the bicritical point accurately. For this 
purpose, we will focus on the region slightly above the 
bicritical point and establish when two separate phase 
transitions are clearly resolved. In this way, we can de- 
termine an upper bound on the bicritical point. 
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1. Signatures of an intermediate paired phase at e — 4.2. 

In order to discern two separate, but close-lying phase 
transitions, we need to establish signatures that can be 
taken as evidence for splitting of a transition line. To 
this end, results are presented for the case when e = 4.2. 
We find unambiguous evidence for two separate phase 
transitions. Remarkably, at smaller system sizes we find 
characteristics of the phase transition consistent with a 
first order transition, and it was interpreted as such in 
Ref. ilH (e = 4.2 corresponds to g ~ 1.88 in the units 
of Ref. [m This Reference gave the estimate for the po- 
sition of the bicritical point at g ^ 2.0.) As we shall 
see, performing computations on larger systems leads to 
a different conclusion. The reason is that finite-size ef- 
fects will disguise the existence of separate transitions 
and make them appear as one. 

In Fig. |4] results are presented for four different ob- 
servables obtained at 12 different system sizes, L e 
{8, . . . , 56}, in a coupling range covering both phase tran- 
sitions. In panel (a), results for the specific heat are 
given. When system sizes are small, it is only possible 
to resolve one peak in the specific heat. However, when 
L = 40, it is possible to resolve a bump to the left of the 
peak. The bump, which corresponds to the 0(3) order- 
ing phase transition, becomes more pronounced when L 
increases. This behavior suggests that there are two tran- 
sitions instead of one. Moreover, in the inset of panel (a) 
we study the scaling of the peak on a log-log scale. When 
L is small, there is a rather steep and slightly increasing 
slope. However, at higher values of L there is a definite 
change in the slope towards smaller values, correspond- 
ing to a sudden slowing down in the growth of the peak. 
This behavior should clearly be associated with resolving 
separate transitions with increasing L. 

In panel (b) of Fig. [4j results for the third moment of 
the action are presented. When system sizes are small, 
it is only possible to resolve a characteristic form cor- 
responding to a single phase transition. However, at 
L > 40, a secondary form is developing to the left of 
the original form, resolving the 0(3) ordering transition. 
When studying the scaling of the quantities (AM3)hcight 
and (AM3)„idth in the insets of the panel, it is clear that 
they both exhibit slope changes associated with resolving 
both transitions.'^ 

The Binder cumulant is given in panel (c) of Fig. |4j 
and its crossings are given in the inset of the panel. 
By considering the crossings with largest L, we find 
that the critical point of the 0(3) ordering transition is 
(3c = 2.347 ±0.001, a value that corresponds well with 
the leftmost transition point in panel (a) and (b). Note 
that there is a non-monotonic behavior in the coupling 
values of the Binder crossings. Hence, by studying small 
systems only, one might be misled to overestimate the 
critical point of the phase transition. 

In panel (d) of Fig. |4j we show results for the quantity 
^PduaiC'lmin): ^^d the Corresponding crossings are given 
in the inset. We estimate the critical point of the Higgs 
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Figure 4. (Color online) Monte Carlo results for four different 
quantities and 12 different system sizes obtained for a coupling 
range covering two separate, but close-lying phase transitions. 
The gauge field coupling e — 4.2. For clarity, the panels 
only show results for L £ {16,24,32,40,48,56}, but insets 
include all 12 sizes, L £ {8, . . . , 56}. Panel (a) shows results 
for the specific heat Cv , and the inset shows the scaling of the 
peak C„,max in a log-log scale. Panel (b) shows the results 
for the third moment of the action M3, and the insets show 
the scahng of (AM3)height and (AM3)width in a log-log scale. 
Panel (c) shows the Binder cumulant U4,, and the inset shows 
the coupling /3cross where the Binder curves cross as a function 
of {LiL2)~^^'^ where Li and L2 are the two actual sizes. Panel 
(d) shows the quantity Lp^l^^i{q^i^) and the inset shows the 
coupling where the curves cross. Lines are guide to the eyes. 
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Figure 5. (Color online) Log-log plot of the value of the spe- 
cific heat peak C„,max as a function of system size L for seven 
different values of eg {3.0,. ..,4.2}. The dotted line corre- 
sponds to the slope expected for a first-order transition, ac- 
cording to Eq. (111. For e > 3.8, the scaling of C„,max shows 



a negative curvature, instead of curving up towards the first- 
order characteristic scaling line. From this, our upper bound 
on the position of the bicritical point in the phase diagram 
would be e = 3.8. Lines are guide to the eyes. 



transition to be /3c = 2.353 ± 0.001 by a crude extrap- 
olation to the thermodynamic limit. Hence, the critical 
point of the Higgs transition is significantly different from 
the critical point of the 0(3) ordering transition. 

The results in Fig. |4] show that it is of particular im- 
portance to simulate large systems in regions where there 
might be multiple phase transitions in multicomponent 
gauge theories. Discarding data points for L > 20, the 
crossings in panel (c) and (d) appear to converge to the 
same coupling. In panel (a) and (b), we would only re- 
solve a single phase transition with rather strong thermal 
signatures. 
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Figure 6. (Color online) Log-log plot of the FSS of the 
height (upper panel) and the width (lower panel) of the 
third moment of the action, for seven different values of 
e £ {3.0,. ..,4.2}. The dotted lines correspond to the slope 
expe cted for a first-order transition, according to Eqs. ( 15 1 
and ( 16 1. Lines are guide to the eyes. 



2. Monte Carlo results for e G {3.0, . . . , 4.6} 

We first turn our attention to the region with e < 4.2 
to look for the signatures that we have established above. 
Fig. [Sjshows the FSS of the peak in the heat capacity for 
e € {3.0, . . . ,4.2}. The results show that there is a def- 
inite change in the slope of the scaling of C^,max, also 
for e — 4.0 and 3.8. Note that this signature of splitting 
appears at higher L when e is reduced, corresponding to 
the coupling difference between the two transitions being 
smaller. The slope of the dotted line in Fig.[5]is the slope 
of a first-order transition [see Eq. (11 1]. For all values of 
e in Fig. [5| we find that for small and intermediate L the 
slope is steep and increasing, and one might be tempted 
to conclude that they all are first-order transitions. How- 
ever, the change towards a smaller slope, that we find for 
large L and e G {3.8, 4.0, 4.2}, is indeed inconsistent with 
a single first-order phase transition. 



In Fig. |6j we show the FSS of (AM3)height and 
(AM3)width- Observe that the same signatures of split- 
ting appears for e S {3.8, 4.0} as found for e — 4.2 above, 
namely that the slope of (AM3)hcight changes to a smaller 
value and the slope of (AM3)width changes to a higher 
value. This is again inconsistent with the scaling of a 
single first-order transition. For a first-order transition 
the slopes should converge towards the scaling for first- 
order transitions, given in Eqs. ( 15 ) and ( 16 ) (see Ref. [T^ 
for an example). 

To determine the positions of the 0(3) ordering tran- 
sition and the Higgs transition, the finite size crossings 
of C/4 and iPduai(*lmin) are given in Fig. Ivlfor eight dif- 
ferent values of e e {3.2, . . . , 4.6}. For e € {4.0, . . . , 4.6}, 
the C/4 crossings and the Lp^^^i{q^ij^) crossings clearly 
extrapolates to difTerent couplings as expected for two 
separate transitions. Also note the corresponding non- 
monotonic behavior for the Binder crossings. When the 
coupling difference between the two phase transitions de- 
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Figure 7. (Color online) P lots of tlie finite size crossings of tiie 
Binder cumulant Ui [Eq. ( 19 1], and the quantity ipduai(lmin) 
for eight different values of e £ {3.2, . . . ,4.6}. The x- values 
are given by (LiI/2)~^^^ where L\ and L2 are the two sizes 
that form the crossing. 



creases, larger systems are needed to resolve this feature. 
For e = 3.8, we observe that the leftmost Ki crossing 
[Li = 80, L2 ~ 96) deviates, consistent with the non- 
monotonic behavior for the larger e values. For the sizes 
available, the crossings seem to converge to the same cou- 
pling value for e £ {3.2, . . . , 3.6}. 

The results in Figs. [Sj |6] and [7j show that there are 
two separate transitions when e > 3.8. We thus estimate 
that the bicritical point must be below e — 3.8. Clearly, 
the system sizes we are able to reach are too small to 
conclusively determine if there are separate transitions 
for e < 3.8. However, in order to estimate the bicriti- 
cal point ebc, in Fig. |8]we show results for the coupling 
difference between the two phase transitions A/3c as a 
function of the coupling e. To estimate when A/3c — > 0, 
in the lower panel, we show A/3c as a function of e — e* 
on a log-log scale where e* is some trial value as labeled 
in the key of the figure. If e* « Cbc, a straight line should 
be expected. A positive curvature suggests that e* > ebc 
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Figure 8. (Color online) Plot of the difference in the critical 
coupling between the Higgs transition and the 0(3) ordering 
transition, A/Jc- A/3c is determined by calculating the differ- 
ence between the I/Pduai(flmin) crossing and the U4, crossing, 
and averaging over four of these differences with largest value 
of (LiI/2)^'^^ (i-e., the four leftmost data points from the pan- 
els in Fig. [7|. We only include results for e > 4.0 where the 
non-monotonic behavior of the Binder crossings can clearly 
be resolved. Upper panel: A/3c as a function of e. Lower 
panel: Log-log plot of A/3c as a function of e — e* where e* 
is given in the key. Positive curvature suggests that e* > ebc, 
negative curvature suggests that e* < etc and a straight line 
suggests that e* « ebc- Lines are guide to the eyes. 



and a negative curvature suggests that e* < ebc- Since 
there is a clear positive curvature both for e = 3.8 and 
e = 3.6, this suggests that Cbc < 3.6. Note that the 
results given in the lower panel of Fig. [H] essentially is 
an extrapolation of the difference A/3c (which also is an 
extrapolation) in the upper panel to find the point ebc 
where A/3c = 0. As it will be clear below, even at the 
largest system sizes accessible for us, we could not prove 
that there is a single first-order transition at e = 3.6. 
Therefore, simulations of even larger systems are needed 
to determine more accurately the existence and the po- 
sition of Cbc- 

Our estimates for the bicritical point differ from the 
results in Refs. [T3] and [TH which studied substantially 
smaller systems. Our upper bound Cbc < 3.8 corresponds 
to /Cbc > 0.151 in Ref. [T31 This means that a part of the 
line that was interpreted as a direct first-order transi- 
tion in that work, in fact are two separate transitions. 
Moreover, the upper bound ebc < 3.8 corresponds to 
.9bc < 1-65 in Ref. [TH where the bicritical point was esti- 
mated io g ~ 2.0. 
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Figure 9. (Color online) Histograms of the probability distri- 
bution of the energy per site H/L^ , for e G {3.4, 3.6, 3.8, 4.0}. 
In every case the flattest (or most bimodal) energy histograms 
were found by reweighting in the vicinity of the pseudocrit- 
ical coupling corresponding to the peak of the specific heat, 
Cu, max. The areas under the curves are normalized to unity. 



3. Signatures of a weak first-order transition 

Although we are led to a different conclusion concern- 
ing the phase diagram than Refs. fT5] and [Til for e > 3.8, 
we find some of the same thermal signatures. As men- 
tioned above (see Figs. [5]and[6]), when systems are too 
small to resolve two phase transitions, the Monte Carlo 
results show that the scaling of C^ max and (AM3)height 
are almost as one would expect for a single first-order 
transition. Moreover, when investigating the energy dis- 
tributions for e € {3.8, 4.0} in Fig.[9j we find that the his- 
tograms are broad. Also, in contrast to previous works, 
we have resolved bimodal structures for e G {3.4,3.6}. 
This could be interpreted as evidence of a first order 
phase transition. At the same time we note that they 
only appear at the largest system sizes. Thus, it is dif- 
ficult to determine if the correct scaling for first-order 
transition is obeyed.'^^'^ The histograms that appear at 
the largest system sizes, have not yet started to evolve 
into distributions resembling delta functions. In particu- 
lar, for the system sizes which we can access the dips in 
between the peaks are still increasing with system size, 
rather than decreasing. The latter is required for drawing 
a firm conclusion that there is a direct first-order phase 
transition at e = 3.4 and e ~ 3.6. Although rare, there 
are examples in the literature where bimodal energy dis- 
tributions are found in cases with no first-order phase 
transition.lSlHlIl 

For e = 3.8, we do not resolve any bimodality, but 
the histograms are wide. The width of the histograms 
decreases and the fiat top structure disappears when L 
increases. This is not consistent with a single first-order 
transition. Note that if this point is located slightly above 
the bicritical point, then according to a mean-fiel d ar- 
gument, the Higgs transition should be first-order !^° l ^^ l 



Also, as mentioned above, the instability of composite 
vortices in type-II SU(2) theory suggests that the system 
should be a type-I superconductor in the proximity of the 
paired phase (since the paired phase results from prolif- 
eration of composite vortices), with a possibility of a first 
order transition via Halperin-Lubensky-Ma mechanism. 
We did not consider large enough system sizes to resolve 
this issue. 

Combining the results in Figs. [5j [6] and |9j it appears 
that for couplings slightly above the estimated bicritical 
point, there are strong thermal signatures in terms of 
broad energy distributions and rapidly increasing peaks 
in the specific heat and the third moment of the action. 
However, when system sizes are larger, we can explicitly 
see signatures of splitting for e > 3.8. We cannot exclude 
the possibility that this may also be the case for some of 
the couplings with e < 3.8. Indeed, the crude extrapo- 
lation in Fig. [8] suggests that e = 3.6 also is above the 
bicritical point. If so, we should expect to see signatures 
of splitting for system sizes larger than those available in 
this work. On the other hand, the strong thermal signa- 
tures we find for e < 3.8 can also be consistent with a 
weak single first-order transition. In that case, we should 
expect to see that proper first-order scaling is obeyed for 
larger system sizes. 

Summarizing this part, we find that the strongest 
signatures for a single first-order phase transition were 
found at e = 3.4 and e = 3.6. Previous works on smaller 
systems did not resolve bimodal structure at these cou- 
plings. For e < 3.4, we did not find any bimodal struc- 
ture in the energy histograms at the system sizes which 
we can reach. 



C. The flowgram method 

To analyze situations where it is difficult to resolve and 
analyze bimodal structures in histograms such as those 
considered above, the authors of Ref. id proposed the 
fiowgram method. By rescaling the linear system size 
L — )• C{g)L, where g = /{4,(3) and where C{g) is a 
monotonous scaling function of the parameter g, it may 
be possible to collapse curves for various physical quan- 
tities computed at the phase transition, for different sys- 
tem sizes and coupling constants, onto a single curve.'23lf 
such uniform scaling is found for all coupling constants, 
one may conclude that a phase transition has the same 
characteristics for all these coupling constants. For in- 
stance, if a first order phase transition were to be found 
for large coupling constants, and the scaled plots fall on 
a single line for all other coupling constants, one may 
conclude that the transition is first order for all these 
coupling constants. To draw such a conclusion, it is very 
important that a broad enough window of systems sizes 
L is considered, such that there is adequate overlap of 
datapoints for all coupling constants, when the data are 
plotted in terms of C{g)L. 

In Fig. 10 we show results of a flowgram analysis 
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of the quantity iPduai(^min) along the 0(3) ordering 
transition Une. For this analysis, the phase transition 
is defined to be at the couphng where the Binder cu- 
mulant U4 — 0.775. With this definition, we will fol- 
low the 0(3) ordering transition line. As mentioned 
above, i/0duai(*3min) is a universal quantity for a con- 
tinuous Higgs transition, whereas it will diverge ~ L for 
a first-order transition. We clearly see such diverging 
behavior when e > 3.6 (not shown here) and the FSS 
is consistent with iPduai(lmin) ^- Refs. [13] and 
1271 this was interpreted as a first-order phase transition. 
However, a diverging Lp^j;^[(q^i,j) is also consistent with 
being above the bicritical point when following the tran- 
sition line of the 0(3) ordering transition. Hence, the 
results in Fig. [TO] correspond well with there being two 
closely separated phase transitions for these values of e, 
see Figs. |4][8] above. 

For e € {3.0,3.2,3.4}, the flowgram analysis suggests 
that iPduai('lmin) diverges, but the FSS is weaker than 
~ L for the sizes available. This is consistent with either 
being above the bicritical point, or with a first-order tran- 
sition. For smaller couplings, the large size behavior of 
the flowgrams is hard to determine. In particular, for 
the couplings e < 2.0 the fiowgrams seem to converge 
slowly to a fixed value, but one cannot rule out diverging 
behavior at larger sizes. 

In Fi g. [Tl] we plot the results for the flowgram data 



in Fig. [101 in terms of the variable C{g)L on a log- 
log scale, using the scaling function C{g) — 3.0324(7 + 
0.0997[exp(4.10055) - l].E2]For large values of C{g)L, the 
collapse appears to be good, and consistent with Ref. [27] 
In our case, we note that for various couplings there are 
sizeable finite-size effects which make it impossible to col- 
lapse smaller systems onto the same master curve. Re- 
moving the data points for the smallest systems for each 
coupling constant would improve the collapse consider- 
ably. 

What can the results of Figs. [TO] and [TT] tell us about 
the character of the phase transition, and about the ex- 
istence of a possible tricritical point separating a line of 
first order phase transitions from a critical line? In Fig. 
[11] the presence of a tricritical point and a line of second 
order phase transition would show up as a bifurcations of 
the master curve at large C{g)L. In Ref. [TO] a tricritical 
point in a global U(l) x U(l) model was detected via a 
breakdown of the curve collapse just below a tricritical 
point. We did not observe such a breakdown of the curve 
collapse for the NCCP^ model. There may exist special 
cases where the universalities of the line of second order 
phase transitions and of a tricritical endpoint are quite 
similar. Then, one may not be able to resolve different 
plateaus at finite system sizes. In such a situation for 
large couplings e, we would have the behavior shown in 
Fig. [11] For small couplings, there should appear an- 
other horizontal branch of the scaling function at large 
values of the argument C{g)L, were a tricritical point to 
exist. The results in Fig. [Tljshow no such feature. How- 
ever, note that the data points for e < 1.50 only extend 
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Figure 10. (Color online) Flowgram of iPduai(Qmin) along 
the 0(3) ordering transition line. In this analysis, the critical 
point is fixed by C/4 = 0.775. Then i/5duai(cimin) measured 
at this point and plotted as a function of system size L. The 
results are given for 13 different values of e G {0.5, . . . , 3.4}. 
The upper panel shows results on a normal scale and the inset 
zooms in on the results for e G {1.0, . . . , 3.0}. The lower panel 
shows the results on a log-log scale. Lines are guide to the 
eyes. 



to about the middle of the plateau in Fig. [TT] This il- 
lustrates the fact, which is also obvious from the lower 
panel of Fig. [TOj that for small couplings e < 1.50, we 
have not reached large enough system sizes to be able 
to ascertain if the curves are horizontal, or if there is 
an upward curvature in any of the curves for e < 1.50. 
Consider for instance the coupling e = 1.50, which is the 
curve in Fig. [TO] which features the most pronounced 
horizontal part for the system sizes we have studied. In 
this curve extends out to C{g)L w 120, which is 
To ascertain whether this 
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Fig. _ ^ 

in the middle of the plateau 
curve falls on the upward curving master curve or contin- 
ues horizontally would require an extension of the curve 
out to C{g)L « 400, or system sizes of about 300'^. No 
computations have been performed on these types of sys- 
tems remotely approaching this range. Another way of in 
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Figure 11. (Color online) Rescaled flowgram of the data in 
The system size L is rescaled by L — >■ C{g)L where 



Fig. 10 



C{g) = 3.0324g + 0.0997[exp(4.1005g) 
Lines are guide to the eyes. 



1] with g = e7(4/3). 
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Appendix A: Mapping the NCCP^ model to a 
J-current model 



principle detecting a tricritical point would be as follows. 
Suppose that one, in order to get good data collapse for 
the entire range of coupling constants would need to re- 
sort to two different types of scaling functions, one below 
some coupling constant and another one above this cou- 
pling constant. At the point where these functions are 
joined, one typically has a non-analyticity. One can thus 
in principal locate a tricritical point at g — gtri by detect- 
ing a non-analyticity in C{g)J^ With our current data we 
have not resolved such a feature in C{g). 



VI. SUMMARY 



In this work, we have studied the three dimensional 
SU(2)-symmetric noncompact CP^ model. We have im- 
plemented an algorithm which permits us to perform an 
investigation of the model at substantially larger system 
sizes than those reached in previous works. It has been 
shown that at couplings e — 3.8 and e = 4.0, which were 
previously estimated to belong to the regime where the 
system undergoes a single first-order phase transition, 
certain signatures should be taken as direct evidence of 
two separate phase transitions. Hence, we conclude that 
a bicritical point must be located below e = 3.8. We find 
bimodality in histograms, consistent with early stages of 
development of a first order transition, at e = 3.4 and 
e = 3.6 (though the histograms do not yet resemble two 
i5-functions and thus indeed it cannot represent a proof 
of a first order phase transitiorP^tSEI! Although our 
estimate for the position of bicritical point is different, 
the data collapse which we find is overall consistent with 
Ref.m 



We start with the lattice formulation of the NCCP^ 
model, 



c,r fi,r 

S = Sf, -\- Sg ^ 

c,/j,,r 

|V'l,r|' + |V'2,r|' = 1 Vr, 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 



where we have introduced t = (3/2 and g = e^/(4/3) - 
the same coupling constants as in Ref. O Writing the 
complex fields on polar form, 



dV'cr dip* 



V'c,r = Per (A6) 

dOc^r / Pc,rdpc,T, (A7) 



we note that the constraint ( A5 ) becomes 

r' = l, Vr, 



(A8) 



which describes the unit circle in the first quadrant of 
the /9i.rP2,r-plaiie (since pc,r > 0). This means that we 
can incorporate the constraint directly into the integral 
by introducing the new field 0, 



Pi^r=COS0r, P2,r^Sinc 



(A9) 



n 



COS 0r sin 0r d(/)r 



(AlO) 
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such that (All, (A3) and (A5) can be replaced by 



It is convenient to introduce the non-negative bond 
subcurrents 



Z = IT / d9i r d92 r / cos (/)r siu (/)r d0r TT / r 6 

(All) 

St = -tY^ [cOS0rCOS0r+yi ( e'(''i-~''i--+'^+^^-) + C.C.) 



-s 



sm (pr sm ( 



-c.c.^ 
(A12) 



Next, we focus on the the ^-dependent part of the in- 
tegrand, namely exp(— S't), aiming at replacing this field 
with a J-current field. First we symmetrize (A12): As- 



suming periodic boundary conditions and using that 



A 



(A13) 



we get 



K.r 

-I- sin0rsin(/)r+ft (e'(''^---"''^ '+~+^-'-) +C.C.) , (AM) 

where k runs over negative as well as positive lattice di- 
rections, K € {±x,±y,zLz}. Then we split exp(— S't) into 
its individual factors and Taylor expand each of them: 



-St 



n E E 

,.=0 /2.'K,r=0 

(I COS (j)r C0S(?!)r+K,) 



fel,».r+'l 



, ■ < \k-2,K,r+l2,K.T 

^2 Sm(/)rSm(/)r+Kj ^ 

=,i(fcl,«,r-il.K,r)(ei,r-01,r + ft+-4,.,r), 



(I 



(A15) 



ai(fe2,»,r-i2,»,r)(e2.r-02,r + ft+A,.,) 



The factors of the product over the lattice and directions 



in ( A15 1 may be rearranged such that all the terms con- 
taining 6'c,r are collected into one, 

Q-St = ^ c''''^,rEJfec,»,r-ic,«.r-fcc,«,r-»-Hic,«,r-s) 

X (Everything else). (A16) 

Here {k,l} denotes the set of all possible Taylor expan- 
sion index field configurations. Inserting this in the par- 
tition function (All), the ^-integrals may now be per- 



formed. The result is Dirac delta functions (up to an 
irrelevant scaling factor, which we ignore) at each lattice 
point, revealing the ("J-current") constraint 

^ ^ ^c,K,r ^c,K,r ^c,K,r— K ^" ^c.K,r— K O5 Vc, T. (Al7) 



Jc,K,r = K.K.v + Ic -K,r+k e No, (AI8) 

as well as the total bond currents 

-^C,K,r — Jc.K,T Jc, — K.X + t^ ^ (-^-^^) 



Reordering the sum, the constraint (A17| then simplifies 
to 



^^c,K,r = 0, Vc,r; 



(A20) 



the current conservation in each component at each lat- 
tice site. 

Getting rid of the 0-field we turn our attention to the 
(/>-field. The terms containing 0r for a given r are on the 
form 



. cos 



sm 



X (Everything else) 



A/'l,r!A/'2,r! 



2(AAi,r+AA2,r + l)! 

where, using ( |A18[ ),( |A19[ ) and ( |A20| ), 



X (Everything else), (A21) 



c,K,r I '^c,K,r— I ''c,K,r— n 



Jc,K,r e No 



K.r-\-K, 



(A22) 



The Taylor expansion (A15) contains an index field 
dependent factor as well. 



En 

{k,l} c,K,r 



~k 1/ ^ 

'^c,K,r •''C,K,r • 



(A23) 



which we want to write as a function of the J-subcurrent 
field instead. It is easy to see that 

by reordering the terms in the product. Using the def- 
inition (AI8), as well as some standard combinatorial 



results, we may rewrite the denominator part of (A23) 
as 



E n 1; 



En E 



1 



k 1; I z-^ 11 I. \( J _ t. V 

=En7^ E (^^^ 



En 

{J} c,K,r 



7 I ' 



(A25) 
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where {J} denotes the set of ah possible subcurrent con- 
figurations. (There is no problem in summing k away, 
as it is an independent variable, and all other terms in 
the partition function are exclusively J-dependent - as 
we will see in a moment.) Inserting ( A24 1 and ( A25 1 into 
(|A23|) gives 



En 



I I' 



(A26) 



which is what we desired. 

Lastly, we want to integrate out the gauge field. The 
gauge field dependent factors of (A15I are on the form 



exp 



c,K:,r 



exp 



c,/i,r 



A T 



(A27) 

Note that the summation is over only positive directions 
on the RHS. (The RHS is found by expanding and re- 
ordering the sum in the exponent on the LHS and ap- 
plying the identity (A13I and the bond current defini- 
tion ( |A19[ ).) Combining ( |A27[ ) with exp(-5'g), the total 



(up to an irrelevant scaling factor) 



Yl / dAr exp^ i 

/i,r /i,r 



OC 



- f E ^c.M.rK.r'/c',M,r'^ , (A28) 



c.c 
/^,r,r' 



where we have applied the Coulomb gauge A^A^ ,. — 0. 
Vr.r' is a long range potential given by by the inverse 
Fourier transform 



E 



sm 



2 fQf^ 



(r - r') , (A29) 



where is the fi component of the Fourier space wave 
vector q. 



gauge field contribution to the partition function reads (A28), leaving out trivial scaling factors, we end up with 



Combining everything, (A20), (A21), (A26), and 



E 



n 



n 



(AAi,r+AA2,r + l)! 



(A30) 



c,c , 



which is a J-current formulation of the NCCP^ model, see also Ref. [HI 
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